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We study the condensation phenomenon for a system of charged bosons in the presence of an 
external magnetic field. We show that condensation happens for a definite critical temperature 
instead of through a diffuse phase transition. The essential ingredient, overlooked in previous anal- 
yses and accounted for in this work, is the treatment of the plasma screening effects by means of 
resummation. We compute the critical temperature, for the case in which the condensate is made of 
charged pions and for typical densities found in compact astrophysical objects, for small and large 
values of the magnetic field. We show that the magnetic field catalyzes the onset of condensation 
at very small and at large values of the magnetic field, and that for intermediate values the critical 
temperature for condensation is lower than for the zero magnetic field case. 

PACS numbers: ll.10.Wx, 67.85.Jk, 26.60.Dd 
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I. INTRODUCTION 



The possibility that a charged pion condensate may 
occur in the interior of neutron stars has been repeatedly 
examined in the past. This possibility is raised by the 
large isospin imbalance between neutrons and protons 
which favors reactions that make neutrons decay into 
negative pions under appropriate conditions. The equi- 
librium thermodynamic conditions obeyed by a pion con- 
densed state in dense neutron matter and in neutron stars 
have been discussed long ago. In particular, Ref. [IJ stud- 
ies the criteria for the appearance of pion condensation 
in neutron matter in terms of the pion Green's function 
(for ageneral review on the physics of neutron stars see 
Ref. Q)- The occurrence of a charged boson condensed 
phase without magnetic fields has also been extensively 
discussed in the literature. Refs. study in-medium 

processes introducing an isospin chemical potential /i/ at 
zero temperature in both phases (\fJ-i\ ^ m w , where ra, 
is the pion mass), analyzing the formation of a charged 
pion condensed phase. This phenomenon was discussed 
in electrically neutral dense quark matter in Refs. @- 
9]. Finite temperature corrections, in the frame of chiral 
perturbation theory, have been considered in Ref. (Ioj |. 
extending the discussion also to other condensates, like 
the chiral condensate or the axial-isospin charge density 
condensate, in Ref. (Tlj . 



The situation becomes even more interesting when con- 
sidering that neutron stars possess large magnetic fields 
whose effects should also be included when studying the 
condensation conditions. Recently, in Ref. [l2j the role 
played by the coupling of ir° to a magnetic field via the 
triangle anomaly has been considered, showing the emer- 



gence of a 7T domain wall for values of the magnetic field 
strength B larger than a certain critical value. This could 
also happen for (77, 77') states when B ~ 10 7 — 10 19 G. 

Magnetic fields can also play an important role in the 
dynamics of systems where charged pions are copiously 
produced, such as relativistic heavy-ion (RHI) collisions. 
Recently, the importance of large magnetic fields for the 
evolution of QCD matter produced in noncentral RHI 
collisions has been discussed in Ref. [lU as well as their 
influence on the phase structure of QCD, with emphasis 
on the chiral symmetry restoration and deconfining tran- 
sitions. In Ref. a discussion of the effective potential 
in the framework of the linear sigma model, coupled to 
quarks and/or Polyakov loop, suggests a richer structure 
of the strong interactions like, for example, a possible 
splitting between chiral symmetry restoration and decon- 
finement in the presence of magnetic fields. The influence 
of the external magnetic field on the formation of CP- 
odd domains in RHI collisions has also been discussed in 
Ref. [l5| . A decrease in the confining critical temperature 
was found in [l6l ]. where a hadron-quark transition was 
studied within the MIT bag model. With the above in- 
gredients put together, the theoretical study of a charged 
boson condensate with a finite chemical potential in the 
presence of magnetic fields becomes even more relevant. 
Although this is an old problem, the results from several 
approaches vary in their conclusions. For instance, it was 
long ago argued that a nonrelativistic Bose-Einstein gas 
of charged particles does not condense in the presence of 
a magnetic field, regardless of how weak the field may 
be . This result motivated the search for conditions 
where condensation could take place with magnetic fields, 
in particular to study whether this could ha ppe n extend- 
ing the number of spatial dimensions [18l42C| . In the 
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nonrelativistic case, treating the dimensionality d of the 
system as a continuous variable, it was shown in Ref. (l8| 
that condensation can happen only for d > 4. For pairs 
of bosons or fermions and in the relativistic case, it was 
shown in Ref. [l9[ that for the case when d is taken as 
an integer, condensation happens for odd d > 5. A sim- 
ilar conclusion was reached in Ref. [2(|, although these 
authors also realized that the lowest Landau level can 
play the role of the ground state to accommodate a large 
charge density in the d = 3 case. 

The common feature of all of the above-mentioned 
analyses is the definition of the condensation condition 
which is taken as the equality of the chemical potential 
and the ground state energy. However, in the presence 
of a magnetic field, this condition leads to a divergence 
of the particle density for that state. Indeed, since for 
a constant magnetic field the energy levels separate into 
transverse and longitudinal (with respect to the magnetic 
field direction) and the former are described in terms of 
discrete energy levels, the divergence of the Bose-Einstein 
distribution when the chemical potential is equal to the 
lowest energy level can only be cured in a larger than 
d = 4 number of spatial dimensions. 

The implications of this condition were recognized in 
Ref. [2l[ where it was argued that when the temperature 
T is much lower than eB one can already consider that 
the system occupies only the lowest Landau level. This 
means that the value for the chemical potential to com- 
pute the ground state density does not need to be equal to 
the lowest energy. In this picture the occupation of this 
state occurs without the need of having a critical tem- 
perature, that is, the system undergoes a diffuse phase 
transition. 

Nevertheless one can argue that if in the absence of 
a magnetic field the system is already in the condensed 
phase with a macroscopic fraction of the population oc- 
cupying the lowest energy level, a slow turning on of the 
magnetic field should not lead to the instantaneous de- 
struction of the condensate. Put in equivalent terms, the 
onset of condensation for small magnetic fields should be 
a phenomenon that takes place at a given critical tem- 
perature T c since it does so in the limit of a vanishing 
magnetic field and the presence of a small one cannot 
drastically change the picture. To implement this idea, 
one should keep in mind that the chemical potential is 
not a number that can arbitrarily be set to take a specific 
value but rather, a function of the thermodynamic vari- 
ables such as temperature and density. Its value should 
be determined by demanding that the ground state is 
populated by a finite charge density. The missing ingre- 
dient that bridges the gap in the analysis is to consider 
the plasma screening effects, which are of course needed 
since we are dealing with infrared phenomena where the 
effective mass is small or may even vanish. 

In this work we study the conditions for the onset of 
a condensed phase for a charged boson system, in the 
presence of an external magnetic field. To mimic the sit- 
uation where there is an isospin imbalance, we introduce 



a finite chemical potential fx. For the description, we re- 
sort to model the boson system in terms of a theory of 
a charged scalar with quartic self-interactions. We show 
that for small and large values of the magnetic field, the 
system presents the magnetic catalysis phenomenon [22j ; 
that is, that the formation of the condensate is favored 
by the presence of the magnetic field. This phenomenon 
has also been found in the context of the Nambu-Jona- 
Lasino model at T = 12311 and in (2+1) dimensions both 
at T = and T ^ [2J], where it was shown that even 
the presence of an arbitrary small magnetic field breaks 
the chiral invariance of the models. A main result of our 
work is to show that when including the plasma screening 
effects, there is a well-defined critical temperature asso- 
ciated with the onset of condensation. A similar calcu- 
lation, using optimized perturbation theory, albeit with- 
out the introduction of a chemical potential, was done in 
Ref. [25[. The authors found that the phase transition 
is always second order and the magnetic catalysis phe- 
nomenon is present for all values of the magnetic field. 
They also found that the critical temperature increases 
with increasing values of the magnetic field. 

The work is organized as follows: In Sec.|TT]we find the 
lowest energy state where condensation happens and de- 
fine the order parameter for the transition. In Scc. lIIII wc 
compute the one-loop corrections to the grand potential 
and set up the discussion for the onset of the condensa- 
tion phenomenon in terms of the existence of a large but 
finite charge density in the ground state. In Sec. IIVI we 
revisit the description of the onset of condensation when 
corrections from interactions are accounted for. We take 
the limit B — > and point out the need to include plasma 
screening effects by means of resummation, even in this 
case. In Sec [V] we explicitly compute the resummed self- 
energy for finite B in the low temperature approximation 
in the limits of small and large magnetic fields. This self- 
energy is then used in Sec. I VII to compute the critical 
temperature for condensation when the charged bosons 
are taken as pions, for typical densities in compact astro- 
physical objects such as neutron stars. We finally sum- 
marize and conclude in Sec. IVIII 



II. ORDER PARAMETER 

We want first to define the order parameter that de- 
scribes the condensation transition. This is a delicate 
task since (as we will show) in the presence of a mag- 
netic field, the condensate docs not correspond to a spa- 
tially uniform state. Let us start by introducing the La- 
grangian representing a charged scalar field <fi with finite 
chemical potential fi interacting with a uniform external 
magnetic field B oriented in the z direction. Working in 
the symmetric gauge, the vector potential corresponding 
to the given magnetic field can be written as 

A = -Bxr. (1) 

2 v ; 
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In Euclidean space and after introducing the magnetic 
field by means of the minimal substitution we get 



= (d T + H)<t>*{0r - H)<j> + |(V - iqA)^ 

+ m 2 \<f>\ 2 + ^ + 5£, 



(2) 



where SC contains the ultraviolet countcrterms and q 
is the (positive) charge associated to the field <t>. Wc 
want to describe the situation where for a given value of 
the chemical potential, the system develops a superfluid 
phase characterized by a boson condensate described by 
a classical held 4> c ; namely, that the field can be expressed 
as 



f 

7T 



(3) 



where <f> is the quantum field, as referred from the real 
classical ground state 4> c . In the presence of an external 
magnetic field, the classical equation of motion does not 
allow a constant value for (j) c [2(| and thus C cannot 
simply be taken as the order parameter for the transition. 
In order to see how one can proceed in such a situation, 
let us find the ground state. Using the Lagrangian in 
Eq. @ we obtain the classical action 

r c = 13 J d 3 J^ c (-V 2 + g 2 A 2 + m 2 - ^ c 



(4) 



where we have discarded a surface term after integration 
by parts and j3 = 1/T. When the space boundary is 
not strictly taken at infinity, the surface term does not 
vanish; however, given the form of the classical solution 
[see Eq. ([7]) below] the boundary contribution can be ne- 
glected for a sufficiently large volume. We first look for 
a solution for the free case (A = 0); thus, the eigenvalue 
problem for the classical equation of motion becomes 

[-V 2 + {qB) 2 {x 2 + y 2 )/4: + m 2 - M 2 ] </> c = £ 2 <p c (5) 

which is recognized as a two-dimensional harmonic oscil- 
lator whose eigenvalues are given by 



£?{p z )=pl + m 2 + {2l+l)qB- l x 2 



(6) 



where I > labels the Landau level. Let us specialize to 
the lowest energy state. This corresponds to I = and 
p z = for which the solution can be written as [27j 



v e 



-qB(x 2 +y 2 )/4 



(7) 



where Vq can be determined from the normalization con- 
dition. The corresponding ground state energy, or effec- 
tive mass squared, is given by 



£o(0) 



qB - n 



(8) 



As anticipated, the solution in Eq. (|7J| is not spatially 
uniform. In order to define an appropriate order param- 
eter, we first normalize the solution over a given spatial 



volume V. This procedure involves finding the average 
over V of d> 2 defined as 



= ~ / 



From Eqs. (J7]) and © one finds 



-*/2*o 



$/2$ 



(9) 



(10) 



where $ = BA is the magnetic flux passing through the 
transverse area A and <j>o = 7t/g is the quantum mag- 
netic flux. Notice that for the analysis, neither A nor B 
can be taken as changing independently but instead that 
<I> should be considered as the relevant variable. Also, 
notice that the term between the parentheses in Eq. (|T0|) 
goes to I as $ — > 0, as expected. The requirement to ob- 
tain the effective mass squared independent of the mag- 
netic flux, leads us to consider <f> c = \/{(j>%) as the order 
parameter to describe the condensation transition. This 
is determined as follows. In terms of <j> c the ground state 
solution reads as 



$/2$ 

I _ g-*/2* c 



1/2 



-qB(x 2 +y 2 )/4 



(11) 



We now look for the value of 4> c that minimizes the classi- 
cal action, this time accounting for the effects of the self- 
interaction (A 7^ 0). Substituting Eq. JTTJ into Eq. (gj) 
we get 



i(< Zj B + m 2 - M 2 )(0 2 ) + A (0 4 ) 



, (12) 



where in general one defines 



1 



-n<E>/4<E> 



n$/4$o 



(13) 



Therefore, the nontrivial minimum of Eq. (|12[) is found 
for a value of the order parameter <p c given explicitly by 



9c0 



/.t 2 — m 2 — qB ( $0 



1 



A V * 

-*/24> \ 1 / 2 



1/2 



1 + e-*/ 2 *" 



qB). (14) 



Notice that for a given value of B (and of $), since 
the thermal occupation of the ground state grows when 
/j 2 — > m 2 + qB, then Eq. (fT4)) means that a macro- 
scopic fraction of the charged particles will occupy the 
condensed state as (j> c o — > + . One can expect that the 
same is true when considering a higher order n in the per- 
turbative expansion of the effective action; that is, that 
the superfluid transition is signaled by the condition that 
4>cn ~ > + . We now proceed to find how this condition is 
realized at order n = 1. 
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III. ONE-LOOP EFFECTIVE POTENTIAL 

The corrections to the value that minimizes the action 
are obtained from the grand potential. For the theory at 
hand, this is given by 



<>( 



D<f>*D(j) e 



(15) 



where S[4>*, <fi] is the action defined in Eq. ([2]) and <j) and 
(j> are related as in Eq. ©. At one-loop order, the grand 
potential has the explicit expression 



^ c ) = (m 2 B - M 2 )(^) + -(^ ) + 



2(3V 



lndct] 



(16) 



where the inverse propagator matrix operator is defined 
as 



-V 2 



- m 

4 



h(h 2 
2 



±s 2 

l ™2 



¥1 



with 



V 2 ± 



(9 r ±/i) 2 + (V ±iqA) 2 . 
Hereafter we use the notation 



ins 



= v™ 2 



qB, 



,(17) 



(18) 



(19) 



for the effective mass in the lowest Landau level. 

The fact that the operators V\ depend on the coor- 
dinates makes it difficult to find the functional determi- 
nant. Nevertheless, since our interest is to explore the 
condensation phenomenon near the phase transition, we 
can resort to expand the grand potential in powers of the 
order parameter near the value <f> c i that minimizes it: 



d 2 n 



The factor 



d 2 n 

Wc 



(20) 



(21) 



can be regarded as an effective mass squared for the 
quantum field <p. Following the discussion at the end of 
Sec. HJ let us explore the behavior of fl((f> c i) and M 2 (4> c i) 
near the supcrfluid phase transition where we expect that 
4>d + . Notice also that M 2 represents the curvature 
of the grand potential in the direction of (f> c . For a given 
value of B (and of <£>) the system is in the normal phase 
when the curvature is positive. As the occupation num- 
ber of the ground state increases, the curvature should 
tend to change sign. Therefore the transition to the su- 
pcrfluid phase is also signaled by the condition M 2 —> + . 
The explicit one-loop expressions for O(0) and Ai 2 (0) are 



n— — oo 



M 2 (0) = m|-M 2 + AT £ j^D 



(22) 



where D is the propagator for a charged scalar in the 
presence of a constant magnetic field. We use the ex- 
pression for D obtained in the Schwinger proper time 
method, given by 



D 



(Is 



cosh(qBs 



\t - \2 i 2 i 2,2 tanh(gBs)l 

e -s[(u„-t / u) +p z +m +p ± g g s \ 



(23) 

where p\ represents the square of the components of p 
transverse to the direction of the magnetic field and u> n is 
a boson Matsubara frequency. Also, in writing Eq. (|2"3")) 
we have ignored a phase factor which does not contribute 
when considering closed loop expressions. Note that the 
one- loop correction in Ai 2 (0) is in fact the self-energy. 
This is not the case when one considers higher loop cor- 
rections. Carrying out the integrations over the trans- 
verse components we get 



M 2 (0) 



where 



qBT 
2tt 



E 

l,n 



dp z 



H + X 



\n(P 2 + to 2 ) 
dp z 



2tt 

qBT x ^ 



2tt 



l.n 



2tt P 2 



P- 



{u n -ill) 2 +p 2 z +qB{2l + l) 



(24) 



(25) 



with I > being the index labeling the Landau levels. 

Let us now separate the T = (purely magnetic field) 
contribution from the thermal dependence by writing 



0(0) 
M 2 {0) 



o 



T,B 



n, 



n 



T.B 



(26) 



where the subscripts T and B denote the temperature 
and magnetic field-dependent contributions. 

Let us first compute the purely magnetic field contribu- 
tions. As shown in the Appendix, these contributions are 
easily obtained using dimensional regularization. They 
can be expressed in terms of the Hurwitz zeta function 
£(s, u) and are given explicitly by 



n B = 



n 



B — 



(4tt) 2 

(*? 

\ TO 

Am 2 
(i^) 2 

(2qB 



1 (2qB\ 

2 ln Uv 



' 2 2qB 



\ TO 



ln<^ r 



2qB) 



/V2tt 



, (27) 



where C'(s, u) = d((s, u)/ds and T(u) is the gamma func- 
tion. In writing Eq. (|27p we have chosen the renormal- 
ization scale in the MS scheme as Aj^g- = me -1 / 2 . With 
this choice, one gets a vanishing contribution in the limit 
B -> 0. 
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As is also shown in the Appendix, the thermal contri- 
bution to f2(0) and M 2 (0) can be expressed as 

qB ^ /•<*> ds e sml -fnVls 



o s 2 1 - e- 2 i Bs 



qB ^ . ,„ , f°° ds e- sm B-P'n*/is 



s 1 - e - 2 9 Bs 



(28) 



Since condensation is a low temperature phenomenon, 
let us approximate Eq. (|28j) in the limit where T <C vtiq . 
In this case, the integrals can be computed using the 
steepest descent method and the result can be expressed 
as 



O 



T,B 



-27rm 4 R T 5/2 



7-^3/2(2) + 



"1 



n 5/2 e n 7 _ I 



n 



T,B 



7Li 1/2 (z) + ^ 



7 J") 



n 3/2 e n 7 _ 2 



where the polylogarithm function is defined as 



Li s (z) 



71 = 1 



(29) 



(30) 



and the fugacity z, scaled temperature r and scaled mag- 
netic field 7 arc defined as 



; (/j-m B )/T 

T 



7 



vtibT 1 



(31) 
(32) 

(33) 



respectively. Notice that the sum over the index n in 
Eq. ([29| corresponds to a sum over Matsubara frequen- 
cies. Also, in writing Eq. (|29]l from Eq. ([28| . we have 
explicitly separated the contribution from the lowest Lan- 
dau level — whose expression is given in terms of the poly- 
logarithm function — from the contribution of the rest of 
the energy levels. This separation proves useful since 
the contribution from the lowest Landau level, unlike 
that from the rest of the levels, is strongly infrared di- 
vergent near the phase transition and must be treated 
separately (2lj |. 

Recall that the charge density is defined as 



(34) 



Near the superfluid transition where A / J 2 (0) — > one can 
compute the charge density by considering only the term 
O(0) in the grand potential, therefore, at one-loop order, 



the thermal part of the charge density is obtained from 
the first line of Eq. as 



m 3 R r 3 / 2 



lLi 1/2 (z) + 



TV) 



T- 

^ n 3 / 2 e™f - 1 



- {p -> -»}■ 



(35) 



Since the phase transition happens when /j ~ ms the 
terms with {fi — > — pi) in Eqs. (|29]l and (|35|) are negli- 
gible, given that they become proportional to powers of 
the factor exp(— (ms+/i)/T) ~ exp(— 2itib/T) <C 1, and 
thus hereafter we ignore them. Notice that with this ap- 
proximation, the second equation of Eq. (|29|) and Eq. (|35|) 
imply that Ht.b an d p are proportional. We have now 
set up the stage to discuss the condensation phenomenon 
in terms of a finite charge density in the ground state. 



IV. BOSE-EINSTEIN CONDENSATION 
REVISITED 

The existence of a critical temperature and a critical 
chemical potential indicates that the system of charged 
bosons reaches a kind of saturation where the occupation 
of the ground state becomes important. This saturation 
leads to the superfluidity phenomenon. From the com- 
putational point of view, the condensation conditions are 
searched for from the values of the parameters that min- 
imize the vacuum energy. Given that the description of 
the onset of condensation may be obscured by the exis- 
tence of infrared divergent quantities and in order to gain 
insight, let us first revisit how these conditions arc found 
at tree and one-loop level. 

At tree level, the condensation condition is given by 
p = ttib- This condition implies the vanishing of the 
mass term M. 2 (0) in the grand potential, as can be seen 
from Eq. (|2l)|) . If we now consider the one- loop correc- 
tion, from the second equation of Eq. (|29[) and from the 
limiting behavior of Li 1 / 2 {z) for /1 — > tub [z ~ 1) and 
fj, < m B (z < 1) 



Lii/ 2 (z) 



^ + C(|) for z<l 
for z < 1, 



(36) 



we see that A^ 2 (0) diverges when (i ^ raj. This behav- 
ior has been interpreted as the impossibility of the exis- 
tence of a superfluid state in the presence of an external 
magnetic field [I?} ■ However, this usual prescription for 
the onset of condensation, is not adequate in the presence 
of an external magnetic field. 

Consider the situation in the absence of a magnetic 
field. Intuitively, once the superfluid phase is established, 
it is difficult to imagine that this can be instantaneously 
destroyed by the turning on of an arbitrary small external 
magnetic field. One would expect that in case the mag- 
netic field destroys the superfluid state, when the field 
is small, the condensed state should be restored for a 
different temperature. Similar considerations were made 
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for the case of a noninteracting gas in Ref. |2l| albeit 
for a high external magnetic field, in the limit T — > 0. 
In such case it was shown that the charge density in 
the normal phase vanishes, i.e., all the charges populate 
the superfluid phase. However, since the physical condi- 
tions require to have a finite charge density, the chemical 
potential — which depends on temperature as well as on 
this charge density — does not reach the value that cor- 
responds to the ground state energy and therefore no 
infrared divergence occurred. It is important to empha- 
size that Ref. [2l| suggests that under such conditions, 
there is no definite critical temperature associated with 
the superfluid transition. 

To continue gaining insight, let us keep on analyzing 
the case with zero external magnetic field. From Eq. (|35|) , 
the charge density for 7 — s- becomes 



/ mT 



3/2 



Li m {e^-™)I T ). 



(37) 



The function Li 3 /2( z ) is well defined for z < 1, although 
its derivative diverges at z = 1. Notice that there is 
no analytical continuation for Li 3 / 2 { z ) which gives a real 
result for z > 1. We have explicitly 



C(|) + V 47r (! - z ) for * % 1 



for z«l. 



(38) 



Thus, the maximum value allowed for the chemical po- 
tential is fiQ = m (z = 1). This result is interpreted 
as the saturation of the boson system that gives rise to 
condensation at this critical value of the chemical poten- 
tial. Since the charge density is a conserved quantity, 
there must be a critical temperature To for which con- 
densation takes place. In other words pq = p{T$, p). The 
condition to obtain the critical temperature comes from 
Eq. (|3"T)) by setting po = m, resulting in 



2tt 

Tn = — 



p 



C( 



2/3 



(39) 



This is the well-known result for the critical temperature 
of a noninteracting boson gas. For temperatures lower 
than the critical temperature, T < To, the gas can be 
separated into two phases: the normal (N) phase and 
the superfluid (S) phase. The charge density splits into 
these two states p = pn + Ps, where the charge density 
in the normal phase pat is defined as the charge density 
p, evaluated at the critical chemical potential pq. In the 
absence of the magnetic field this is given by 



fmT\ 3/2 , , . 
PN={-^) C(3/2). 



(40) 



The superfluid charge density corresponds to the differ- 
ence between the total charge density and the charge den- 
sity in the normal phase, ps = p — Pn, for a fixed total 
charge density p. 

When interactions are accounted for, the situation 
changes. The condition for the phase transition is once 



again looked for from the vanishing of the effective mass 
squared Ai 2 (Q). From Eq. and the second equation 
of Eq. (|29|) . in the limit B — > 0, the effective mass squared 
is given by 



M 2 (0) 



m — p 
Am 2 / T 



3/2 



2 \2TTmJ 



(H- m )/T\ 



Xm? 



T 

2-7TTO 



3/2 



C(3/2), 



(41) 



which shows that M 2 (0) cannot vanish for p < in. More- 
over, since the function Li 3 / 2 ( z ) cannot be analytically 
continued to real values for z > 1, there is no physical 
solution that sets Ai 2 (Q) = even if we were to consider 
p > m. 

The above results show the need of an extra ingre- 
dient already for B = 0. In this case, the solution is 
well known: since the physical conditions require the ef- 
fective mass squared to vanish, plasma screening effects 
need to be accounted for by means of resummation. Ex- 
amples of the importance to include resummation effects 
have been recently discussed for systems subject to the 
influence of an external magnetic field. For instance in 
Ref. [28| the temperature dependent effective potential 
for a scalar theory, similar to the case here discussed, 
was considered. A resummation of ring diagrams turns 
out to be extremely important to understand the appear- 
ance of a first order phase transition. The scenario has 
also been considered in other theories such as the lin- 
ear sigma model [29l | and the standard model, during 
the electroweak phase transition [30j . In all these cases, 
a second order phase transition turns into a first order 
one. 

In the present context inclusion of resummation ef- 
fects means that the self-energy should be computed self- 
consistently as 



n 



n = XT E / 7S3 D ^ 2 ™ 2 + n)- ( 42 ) 



To simplify the calculation, and in case the coupling A 
is not too large, it is customary to substitute inside the 
argument of the propagator in Eq. (|4"2")l , ft — > Hi. For the 
low temperature expansion for B = and when p = to, 
we see from the second equation of Eq. (|29| and Eq. (|38|) 
that in this case 



(43) 



When the coupling A is not small, a full self-consistent 
treatment is needed. 

In the following section we explain in detail the re- 
summation procedure for a finite external magnetic field. 
From there, the case treated in this section is obtained 
as the limit with B —> 0. 
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V. RESUMMATION AT FINITE B 

As was mentioned in the previous section, in order to 
consistently compute the critical temperature when inter- 
actions are accounted for, it is necessary to consider re- 
summation. Inclusion of these effects allows us to find the 
chemical potential, beyond its tree level value (i = m, by 
consistently accounting for plasma screening. The need 
for resummation is even more dramatic for B ^ 0, since 
when ignoring screening and /i — > mj the charge density 
and the thermal contribution to the self-energy both di- 
verge in the infrared. In Ref. [3l[ the thermal effective 
potential in the A</> 4 theory was computed by a resum- 
mation of the ring diagrams. Effectively, the resumma- 
tion was done through a solution of the renormalization 
group equation. Resummation results in a shift of m 2 in 
the one-loop correction of the grand potential 

fi(o) -> n r (o) = n B + o T ,s| m 2^ m 2 + n 

M 2 (0) -^M 2 ,(0) =m% -n 2 

+ n B + n T , B | m2 ^ ro2+fl , (44) 

where the subscript r is to emphasize that resummation 
effects are included. As a result, the charge density be- 
comes 



dQ r {0) 
9/i 



(45) 



The resummed self-energy II is obtained self-consistently 



n = a 



qBT 
2tt 



E 



l . n 

{T^O}, 



dp z 



2tt P 2 



n 



(46) 



where P 2 is defined in Eq. (|2"5"j) . We consider only 
the thermal contribution, as indicated by subtracting the 
{T — > 0} term. Also since we are interested in computing 
this self-energy near the phase transition, we consider its 
value for (i = itib ■ The self-consistent equation is reduced 
to finding the solution to the expression 



n = n 



T,B\ 



(47) 



For the rest of this section, we concentrate on finding an 
explicit expression for II in some limits. 

In order to simplify the self-consistent equation, recall 
that since we want to include the effect of thermal fluctu- 
ations coming from the resummation of the so-called ring 
diagrams for small temperatures, we use the low temper- 
ature approximation for Ht,b given in Eq. (I29|) . We can 
consider II <C m 2 up to the leading order in the right- 
hand side of Eq. (14"T1) . At leading order, we can set n = 
in all the terms, except in the argument 



; (, I -V m2 +<? B )/' r 



■> 2 ->m 2 +n, 



(48) 



of the polylogarithm function L 1 / 2 (z), since this diverges 
when fi = n%B in the limit II — > 0. Therefore we set 



s ("IB - \/m 2 +n+qB^ /T 
D -n/2m B T 



(49) 



We explore the situation where II <C rriBT, which is 
a good enough condition to control the infrared diver- 
gences. Then we can expand the fugacity as 



1 - U/2m B T. 



(50) 



Using the expression for Li 1 / 2 in Eq. (|3l))) we obtain a 
simplified expression for the self-consistent equation 



II « ^mlr 



3/2 



IkwibT 



1 



n 



n 3/2 e n 7 _ I 



(51) 



This is a third order equation in the variable vll. It 
can be solved explicitly, though the solution is given in 
terms of complicated expressions. It is therefore more 
instructive to find the solution for different values of T 
and B, in particular, the next-to-leading order correction 
in A for the case qB <C itibT here treated. For this 
purpose, we expand Eq. (|5"Tj) for 7—^0. Also, since we 
are close to B = we can write IIi, given in Eq. (|43[) . 
instead of II in the right-hand side of Eq. (jBTj) , resulting 
in 



n 



\ x l 2 qB 



( T ^ 


1/4 Am 2 


( T \ 


\ 2irm J 


+ 2 


\2nm 



3/2 



(52) 

It can be shown that the first term on the right-hand 
side of Eq. (|52[) is the leading correction coming from the 
lowest Landau level. Notice that this term also makes 
explicit the breaking of the perturbative regime as it is 
proportional to A 1 ' 2 . 

We emphasize that resummation is needed only for the 
case when considering small values of the external mag- 
netic field, compared with temperature. In fact, notice 
that the nonthermal contribution to the self-energy TIb 
in Eq. (|27p is negative. Its limiting values can be approx- 
imated as 



XqB 
' (4tt) 2 



In 2 for qB » 77T 



(53) 



Therefore, in this case there is always a sufficiently small 
temperature such that lis + ^-t,b < for \x < m# and 
the effective mass squared A'f(O) 2 in Eq. (|26p can thus 
vanish. 



VI. CRITICAL TEMPERATURE 

Having solved the infrared divergence problem by the 
introduction of the resummation procedure in the self- 
energy, we now proceed to find the critical tempera- 
ture for the supcrfluid phase transition and explore how 
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this critical temperature is modified as a function of the 
strength of the magnetic field, compared to the zero field 
case. As we pointed out in Sec. HV) the charge density 
is a conserved quantity and therefore it cannot diverge. 
From the resummed version of Eq. ([35]), that is, with 
the mass shifted as m 2 — > in 2 + II, we can express the 
chemical potential in terms of the charge density and the 
other parameters, namely, p = p(p,T, B). For a fixed 
charge density and a given value of the magnetic field, 
the critical temperature T c (p, B) can be found from the 
condition that sets the grand potential up to its minimal 
value. Recall that for B = it is well known that the 
phase transition is of the second order. Even more, in the 
case of zero chemical potential, the transition is always 
of second order (25|. Here we will assume that, for B^O 
the phase transition continues being of the second order, 
which is the expected behavior if we slowly turn on the 
magnetic field. 



1.05 



M 2 r (0)\ T=Tc 
dfl r {0) 



dp 



= 

= P, 



(54) 



where p c (p, B) = p(T c , p, B). We proceed to numerically 
solve Eq. ([541 . 

In the low temperature approximation that we are 
considering, the thermal contribution to the self-energy 
is proportional to the charge density, as can be seen 
from Eqs. ([29)) and ([35]) . upon neglecting the terms with 
{p — > —p}. Therefore we can write 



n 



T,B 



2m B 



(55) 



The same is true after including rcsummation effects. 
Since p is constant, we can use the above result to imme- 
diately solve for the critical chemical potential from the 
first line of Eq. (|5"4")) , resulting in 



n 



1/2 



(56) 



Notice that now, finding the critical temperature be- 
comes equivalent to solving Eq. (|35|) with the replacement 
m 2 — > m 2 + fl, where for T and p use is made of their 
critical values T = T Cl p = p Cl with p c given explicitly in 
Eq. (f56|) . Let us now find T c numerically. For this pur- 
pose, we consider a charged pion condensate in the core 
of compact stars. The mass is given by m = m T 140 
McV. The typical surface temperature for a cold neutron 
star is T ~ 100 eV. Since the critical temperature T c for 
small B is expected to be close to the value T c q at zero 
magnetic field, we can use Eq. (f3T)]) to obtain the charge 
density that is needed for the onset of condensation at 
these temperatures. 

Figures [T]and [2] show the critical temperature T c scaled 
to the critical temperature at zero magnetic field T c q as a 



0.95 
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FIG. 1. Critical temperature T c scaled to the critical temper- 
ature T c o as a function of the magnetic field strength scaled 
to m|. The calculation is valid for the case when qB <C itlbT. 
The upper panel corresponds to a density p = 10~ 10 m L l 
whereas the lower panel to a density p = lQ~ 7 m%. 



function of the magnetic field strength scaled to m 2 ^. In 
both cases we consider two different values of the charge 
density, p = 10 _10 m^ and p = 10~ 7 m^, which corre- 
spond roughly to temperatures To ~ 100 eV and To ~ 10 
kcV, respectively. 

Figure [T] corresponds to the case qB <C m^T. In prac- 
tice we consider that the parameter 7 defined in Eq. (|33|) 
is such that 7 < 10~ 3 . Notice the smooth rise of the 
critical temperature as the magnetic field increases. This 
signals that when the magnetic field is slowly turned on 
it catalyzes the formation of the condensate; that is, the 
critical temperature is larger than for the B = case. 
Nevertheless this behavior disappears for higher values of 
the magnetic field and the formation of the condensate 
appears at lower temperatures. We also mention that 
catalysis in this case can occur for a larger range of mag- 
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FIG. 2. Critical temperature T c scaled to the critical temper- 
ature T c o as a function of the magnetic field strength scaled 
to m%. The calculation is valid for the case when qB 3> uibT. 
The upper panel corresponds to a density p = 10~ 10 m^. 
whereas the lower panel to a density p = 10~ 7 m^. 



netic field strengths when the charge density is larger, as 
well as for larger values of the coupling constant A. Fig- 
ure corresponds to the case qB ^> vtlbT . We explicitly 
consider 7 > 10 3 . For these larger values of B the critical 
temperature grows. This means that at some intermedi- 
ate values of qB between the large and small regions here 
considered, the magnetic field catalyzes again the forma- 
tion of the condensate. The catalysis is enhanced as the 
coupling constant A increases, and contrary to the low 
magnetic field case, it is also enhanced for lower densi- 
ties. 



VII. SUMMARY AND CONCLUSIONS 



In this work we have studied the condensation phe- 
nomenon for a system of charged bosons in the presence 
of an external magnetic field. Contrary to what is com- 
monly believed, we have shown that condensation hap- 
pens at a definite critical temperature. The missing in- 
gredient, overlooked in previous analysis and accounted 
for in this work, is the treatment of the plasma screen- 
ing effects by means of resummation. We have explicitly 
computed the critical temperature, for typical densities 
found in compact astrophysical objects, for small and 
large values of the magnetic fields. We have shown that 
the magnetic field catalyzes the onset of the condensa- 
tion at very small and large values of the magnetic field, 
agreeing in these regions with the case of zero chemical 
potential [25|. For intermediate values of the magnetic 
field, the critical temperature for condensation turns out 
to be lower than in the B = case. 

Recall that although the term magnetic catalysis usu- 
ally refers to an enhancement of dynamical symmetry 
breaking by an external magnetic field, the phenomenon 
seems to be universal, appearing in different physical sce- 
narios [32l | such as the one treated in this work, namely, 
a condensate of scalar particles, as opposed to the more 
usually studied case of a fermion condensate. 

We should also mention that the problem has been 
studied by lattice methods as well. The first lattice simu- 
lation for dcconfincmcnt and chiral symmetry restoration 
for two-flavor QCD, in the presence of a magnetic back- 
ground, was done in Refs. |33l. |34|. The result was that 
the transition temperature significantly decreased with 
increasing magnetic field. A similar conclusion was pre- 
sented in [3.5]. In that work it was found that in the 
chirally broken phase, the chiral condensate increased 
monotonically with a growing magnetic field strength. 
In fact, in the chiral limit this behavior started linearly. 
In the same limit and in the chirally restored phase, the 
condensate vanished independent of the strength of the 
magnetic field. On the other hand, in Refs. [33, HI| the 
effect of an external magnetic field on the finite tempera- 
ture transition of QCD was considered. Thermodynamic 
observables including the chiral condensate and the sus- 
ceptibility were measured. The result was that the tran- 
sition temperature significantly decreased with increasing 
magnetic field. Such discrepancies can be originated by 
the fact that in Refs. I33I435I the pion had a large mass. 
In the case of Refs. [37L l38l | light pions as well as an 
improved lattice were used. These seemingly contrast- 
ing results call for a closer look at the phenomenon, in 
particular for the case of intermediate values of the mag- 
netic field strength. This is work that we are currently 
pursuing and will be reported elsewhere. 
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Appendix 

The grand potential for a scalar field in the presence 
of a uniform magnetic field evaluated at <p c = 0, is given 
by 

= ^ trln (^+ +«i 2 ) + itrln(X>_ +m 2 ), 



= fdm 2 T Y. J <g^, 



(A.l) 



n— — oo 



where the operators T>± arc defined in Eq. (JTSJ) and the 
propagator D expressed in the Schwinger proper time 
formalism is defined in Eq. (f2"3"]l . The validity of the 
second line in Eq. (|A.1[) is more easily seen from the 
representation of the grand potential in terms of Landau 
levels, Eq. (|24l) . Thus, to obtain the grand potential, 
we need to compute the term proportional to the one- 
loop self-energy and then integrate over m 2 , neglecting 
unphysical ultraviolet divergent constants. 

The sum over Matsubara frequencies uj n = 2irnT gives 
rise to the Jacobi thcta function 



-TTXTl +2-KZTI 



n— — oo 



which obeys the inversion property [3j 

TTZ 2 /X / 1 

6 a ( z ,x) = =-=- e ' 



Identifying z = 2iTfis and x = 4irT 2 s allows us to rewrite 
the full propagator as 

D _j_ r d S e- s ^ 2 +p 2 z+pi- 



\Fk Jo V 7 * cosh(e_Bs) 

1 °° t) 2 ™ 2 
- + e ~ cosh(/3/m) 



(A.2) 



The ultraviolet divergent term is the first term inside the 
square bracket. The one-loop self-energy, including the 



counterterm 8m 2 , is given by 



n — — rv-i u v / 



n— — OO 



n B + n TiB , 



(A.3) 



where we separate the thermal contribution Ht.b from 
the purely magnetic field contribution lis- This last 
contribution contains the ultraviolet divergence, which 
is canceled with the counterterm 8m 2 using dimensional 
regularization 25]. Introducing the scale factor A, the 
purely magnetic contribution is given by 



n, 



8m 



X 



d d P K s-d 



4^ J (2tt)' 



ds 



r 2 2 i 2 tanh(gBs) i 

-s[m +p z +p ± g s ] 



cosh(qBs) 



(A.4) 



Integrating the p± and p z momentum components and 
expanding in powers of e~ qBs , the purely magnetic con- 
tribution to the self-energy is expressed in terms of a sum 
over the Landau levels 



n, 



8m 



= 8m 2 



2XqBK 



d-3 00 



ds- 



-s(m 2 +qB(2l+l)) 
s (d-l)/2 



4ttA^ 



2XqB ( 
f (An) 2 \ 2qB 
3-d 1 to 2 



2 x (3-d)/2 



2qB 



3-d 



(A.5) 



where C, is the Hurwitz zeta function 

oo 



(A.6) 



Expanding in the number of dimensions 3 — d we get 



n, 



5m 2 



Am 2 



(4tt) 2 

2eB • f 1 m 2 \ 
— C K + 



nr 



3-d 

2 

2eB ) 



7 + In 



4ttA 2 
2qB 



(A.7) 



where £'(s,u) = d s ((s,u). The divergent term is re- 
moved using MS scheme. In order to have a vanishing 
contribution to the self-energy in the absence of thermal 
and magnetic effect, one can choose for the scale factor 
the convenient value A = me -1 / 2 . Using the relation 
£'(0,u) = ln(r(u)/v / 27r) we finally arrive at the expres- 
sion for U B in Eq. (|2"7|) . 

The expression for £Ib in Eq. (|27p is obtained by in- 
tegrating TLb/X with respect to m 2 in Eq. (|A.7[) . and 
using the same value for the scale factor to remove the 
divergent terms in the MS scheme. 

The finite temperature part is given by the second term 
inside the square bracket of Eq. (|A.2[) . After integration 
over p± and p z , the function Ht,b can be written as 
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in Eq. (|28|) . Consequently, by integrating IIt,_b/A with 
respect to m 2 we get the function £It.b hi Eq. 



In the limit T<ms, one can use the steepest descent 
approximation [39|. By scaling the variable s — > s 1 /tiibT 



in Eq. (|28|) . both integrals can be expressed in terms of 

/2^3(so)e /3mB/(so) 



ds'e pmB}(s ' ] g{s') 



\pm B f"(s )\ 



1/2 



(A. 



with f(s) = — (s + n 2 /As) and with the saddle point 
,s = n/2. With this approximation we arrive at Eqs. (p?5)) 
and (J2SJ). 
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